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Abstract. Consider a homogenous fluid membrane described by the Helfrich- 
Canham energy, quadratic in the mean curvature of the membrane surface. The shape 
equation that determines equilibrium configurations is fourth order in derivatives and 
cubic in the mean curvature. We introduce a Hamiltonian formulation of this equation 
which dismantles it into a set of coupled first order equations. This involves interpreting 
the Helfrich-Canham energy as an action; equilibrium surfaces are generated by 
the evolution of space curves. Two features complicate the implementation of a 
Hamiltonian framework: (i) The action involves second derivatives. This requires 
treating the velocity as a phase space variable and the introduction of its conjugate 
momentum. The canonical Hamiltonian is constructed on this phase space, (ii) The 
action possesses a local symmetry - reparametrization invariance. The two labels 
we use to parametrize points on the surface are themselves physically irrelevant. This 
symmetry implies primary constraints, one for each label, that need to be implemented 
within the Hamiltonian. The two lagrange multipliers associated with these constraints 
are identified as the components of the acceleration tangential to the surface. The 
conservation of the primary constraints imply two secondary constraints, fixing the 
tangential components of the momentum conjugate to the position. Hamilton's 
equations are derived and the appropriate initial conditions on the phase space variables 
are identified. Finally, it is shown how the shape equation can be reconstructed from 
these equations. 
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1. Introduction 

A biological membrane consists, in a first approximation, of a bilayer of phospholipid 
molecules in water QjJ |2] . The thickness of the membrane is determined by the length 
of the hydrocarbon tails ~ 50 nm, whereas its overall size can reach 100 //m - it is thus 
sensible to describe it on mesoscopic scales as a two-dimensional surface. Remarkably, 
many of the physical properties of the membrane are completely described by the 
geometrical degrees of freedom associated with this surface. The language that is most 
appropriate to describe the membrane is differential geometry. 

The phospholipid molecules move freely within the membrane so that the membrane 
behaves as a two-dimensional fluid. In a geometrical model, this feature is implemented 
by associating an energy, invariant under surface reparametrizations, with each possible 
shape of the membrane. One should think of reparametrization in its active sense, as the 
freedom of a point to move about on the surface, rather than as a change of coordinates. 

In the early seventies, it was recognized that the leading term in the energy is 
proportional to the cost of bending the membrane - quadratic in the mean curvature 
IU IH]- This geometric invariant is known in the mathematical literature as the 
Willmore functional [6 . The geometry of a fluid membrane will generally also be subject 
to constraints: its area and (if closed) its volume are constant. Moreover, the bilayer 
structure is captured by the area difference between the layers, proportional to the 
integrated mean curvature. Its value is constrained in the minimal curvature model for 
fluid membranes [Zj. This is the model we will focus our attention on. Other curvature 
models have superseded this simple model (see e.g. [S] for a review). In particular, it 
is worth mentioning the ADE model, with a non-local constraint involving the square 
of the area difference [011101 \H\ Our considerations can be extended to this and other 
geometrical models. 

We are interested, in particular, in the equilibrium shapes of a closed fluid vesicle 
minimizing the Helfrich-Canham energy, defined below by Eq.@, which describes the 
minimum curvature model. Typically, the addition of constraints, will increase the 
number of equilibrium solutions. The vanishing of the first variation of the energy 
is expresssed in a single 'shape' equation O H31 El- It is a fourth order non- 
linear PDE. Beginning in the midseventies, with the work of Deuling and Helfrich on 
axisymmetric configurations |15| . considerable effort and ingenuity has been dedicated 
to characterizing these shapes. For axisymmetric configurations, it is relatively 
straightforward to solve the shape equation: the symmetry reduces it to a 3rd order non- 
linear ODE. Even in this simple case, however, solving this equation involves numerical 
approximations [IB- . For more general configurations, analytical techniques have been as 
good as useless outside of perturbation theory. Rather than attempt to solve the shape 
equation itself, one has relied on a combination of computational techniques, among 
them Monte-Carlo simulations and dynamical triangulations to minimize the energy 
(see e.g. PSHEj). 

In principle, how does one go about solving the shape equation when axisymmetry 
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is relaxed? The method used to solve the axisymmetric shape equation numerically 
suggests a strategy. What one does is decompose the higher order shape equation 
into a set of coupled first order equations. Initial conditions for these equations are 
specified, that will need to be tuned to obtain a geometry which closes smoothly. In 
practice, modulo the step involving tuning, what one does is formulate the problem 
in Hamiltonian terms. In this paper, we provide a novel formulation of the shape 
equation, in all its glory, as a Hamiltonian sytem. The basic idea is borrowed from the 
Hamiltonian formulation of General Relativity, and more specifically, from the recent 
Hamiltonian formulation of the dynamics of relativistic extended objects jTHj. Take any 
curve on the membrane surface. The surface is then generated by the evolution of the 
curve along some direction with appropriate initial data. The evolution is determined 
by a Hamiltonian that is derived from the Helfrich-Canham energy, reinterpreted as 
an action for the curve. To do this in one step is daunting technically. The reader 
may be interested in first working through the details in a simplified setting. The 
axisymmetric case has been considered in ^H], where the shape equation is formulated 
as a Hamiltonian system describing the motion of a particle, rather than a curve, in two 
dimensions. While the canonical structure itself is illustrated clearly in this restricted 
setting, it is too simple to capture the geometry. 

The paper is organized as follows. In Sect. 2, we introduce the geometry for 
the surface that describes the fluid vesicle, and the Helfrich-Canham energy. Next, 
we consider an arbitrary parametrized closed curve on the surface, and decompose the 
surface geometry along a tangential basis adapted to the curve: one vector is tangent 
to the curve, the other, not necessarily orthogonal to it, is parameterized by a second 
parameter t. The latter tangent vector acquires an interpretation as the velocity of the 
curve, parameterized by the time t. By implementing this decomposition of the surface 
geometry within the Helfrich-Canham energy, it can be written as an action describing 
the evolution of the curve. The corresponding Lagrangian depends not only on the 
position of the curve and its velocity, but also on its second derivative with respect 
to the parameter t: its acceleration. A second key feature is the presence of a local 
symmetry of the action, reparametrization invariance. Before looking at lipid vesicles, 
in Sect. 3 we consider the simpler, but still non-trivial, case of a soap bubble. The energy 
is proportional to the area of the bubble with a constraint on the enclosed volume. The 
action does not involve second derivatives and we can address the problems associated 
with reparametrization invariance without the added complication of higher derivatives. 
In Sect. 4 we return to fluid vesicles described by the Helfrich-Canham energy, and 
we identify the phase space variables which are appropriate when the Lagrangian 
involves the acceleration; it becomes necessary to introduce canonical momenta not 
only conjugate to the position of the curve, but also to its velocity. We show how these 
momenta can be identified by examining the response of the energy to deformations of 
this curve. The momentum conjugate to position will be shown to be related directly to 
the conserved stress introduced by two of us in [T3j. Not surprisingly, the stress plays a 
role in the Hamiltonian formulation; stresses get transmitted across closed curves j2Hj- in 
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Sect. 5, we construct the canonical Hamiltonian in terms of a Legendre transformation 
with respect to both position and velocity. This Hamiltonian will generate the evolution 
of the curve in terms of the phase space variables. At a Hamiltonian level, the existence 
of a local symmetry, in our case reparametrization invariance, manifests itself in the 
appearance of constraints on the phase space variables at fixed values of the evolution 
parameter. We find two primary constraints, that need to be added to the canonical 
Hamiltonian to generate evolution, as well as two secondary constraints. These four 
constraints are conserved in time. In particular, one of the secondary constraints is the 
vanishing of the canonical Hamiltonian - the hallmark of reparametrization invariance. 
The secondary constraints are the generators of reparametrizations. In Sect. 6, we 
derive the Hamilton equations for this system. Appropriate initial conditions on the 
phase space variables are discussed. Reparametrization invariance manifests itself in 
the appearance of two arbitrary functions in these equations. At the end of the day 
the shape equation has been dismantled into its Hamiltonian building blocks: four 
constraints at fixed value of time and four first order, in time, evolution equations. In 
Sect. 7, we show how the shape equation can be reconstructed using these blocks. We 
conclude in Sect. 8 with some final remarks. Various technical details are collected in a 
set of Appendices. 



2. Helfrich-Canham energy 

We model a fluid membrane as a two-dimensional surface E, defined locally by the 
embedding 

x = X(C) , (1) 

where x = x l are local coordinate in space, £ a local coordinates on the surface E 
... = 1, 2, 3; a, b, . . . = 1, 2), and X = X 1 the embedding, or shape, functions. We 
denote by e a = <9X/<9£ a the two tangent vectors to the surface. The metric induced on E 
is given by their inner product, g ab = e a -eb, and its inverse is g ab ; indices are lowered and 
raised with g ab and g ab , respectively. The area of the surface E is A = J dA = J d 2 !;y/g, 
with g = detg a b- The unit normal n to E is defined implicitly by e a ■ n = 0, n 2 = I. The 
extrinsic curvature tensor is K a b = — n ■ <9 a <9&X, and the mean curvature is K = g ab K a b. 
In terms of the principal curvatures, {01,02}, we have K — c\ + c^- The intrinsic 
scalar curvature can be given in terms of the extrinsic curvature via the Gauss-Codazzi 
equation as 1Z = K 2 — K ab K ab ; it is twice the Gaussian curvature, so that in terms of 
the principal curvatures TZ = 1c\C2- 

The shape of a fluid membrane is obtained by minimizing the Helfrich-Canham 
energy, quadratic in the extrinsic curvature. In the bilayer coupling model, the 
membrane is subject to various global constraints: the enclosed volume, the 
well as the integrated mean curvature (the area difference between the layers) are fixed. 
One then searches for stationary points of the Helfrich-Canham energy F[X] defined by 

F[X] = ^JdAK 2 + f3jdAK + aA-PV. (2) 
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The constant k is the bending rigidity; /3, a and P are Lagrange multipliers enforcing the 
constraints of constant integrated mean curvature , constant area and constant enclosed 
volume V, respectively. The constant a is the bare surface tension, and P is the osmotic 
pressure. Of course, the inclusion of the volume term requires that the surface be closed. 
The volume can be written surface integral: 

V = - [ dA n • X . (3) 
3 J 

Despite appearances, this expression for V is translationally invariant. This is because 
/ dA n ■ a = for any constant vector a on any closed surface. We have not included in 
the energy a term corresponding to the Gaussian bending, i*b[X] = kq j dATZ, with 
the Gaussian bending rigidity, since for a two-dimensional surface, by the Gauss-Bonnet 
theorem, it is a topological invariant, and so it does not contribute to the determination 
of equilibrium configurations. 

The vanishing of the first variation of the energy (j2j), with respect to infinitesimal 
variations of the shape functions X — > X + <5X, gives the shape equation [121 HHJ EH I2T] 



_V 2 iT -y(^ 2 -211) 



+ (3K + o-K - P = , (4) 



where V 2 = g ab V aVb denotes the surface Laplacian, and V a the surface covariant 
derivative. This fourth order non-linear PDE determines the equilibrium configurations 
of fluid vesicles. Note that there is only one equilibrium condition. Naively varying 
X, one would have expected three, since there are three independent deformations. 
However, reparametrization invariance tells us that two of the equations, corresponding 
to tangential deformations, must vanish identically |22j- The only physical deformations 
that enter in the determination of equilibrium are those normal to the surface. 

Our basic proposal in this paper is to consider the surface £ as the evolution of 
some closed curve C. If the curve is parametrized by u, and its evolution by t, the 
surface generated by C will be given by x = X(tt, t). For any fixed t, the tangent vector 
to the curve is X' = <9 U X = dX/du. The infinitesimal arclength along the curve is 
ds = h l l 2 du, where 

h = X' ■ X' (5) 

is the one-dimensional metric along C; the unit tangent vector is t = /i _1 / 2 X'. The 
curve C evolves along the surface vector field X = <9 t X = dX/dt. This vector will not 
generally be orthogonal to C. Let 1 denote the unit vector normal to the curve on the 
surface 1 • X' = 0. We expand X with respect to the vectors 1 and X': 

X = M + JV||X'; (6) 

in analogy with the terminology used in general relativity, we will refer to the 
components N, Nu with respect to this basis as the lapse and shift functions, respectively. 
In terms of X' and X, they are given by 

iV 2 = X 2 — (X " Xy , iV || = ^ X -^. (7) 
h h 
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We now decompose the geometry of the surface X along the basis of tangent vectors 
adapted to the evolution of the curve, {X, X'}. For the induced metric we have 

(N 2 + hN 2 hN\\\ 
9ab = { hN, h ) ; (8) 

and for its inverse it follows that 

9 ab = N~ 2 ( \ T ~f" AT2 ) . (9) 

y \-N\\ N 2 ^ 1 + N 2 J w 

Note that the component g uu can be written also as g uu = X 2 /h. The determinant of 
the induced metric g a t takes the simple form 

g = N 2 h. (10) 

The vectors t, 1, together with the normal to the surface n form a spatial basis adapted 
to the curve. We choose an orientation so that n = t x 1. In terms of X' and X, 

n=^^X'xX. (11) 

The decomposition of the extrinsic curvature along the basis {X, X'} is given by 
/n-X n-X' \ 

** = -(».* n-X")' < 12 > 

which depends on the acceleration of the curve X projected along the normal to the 
surface. For the mean curvature it follows that 

K = g ab K ab = ^(-n-X + j) , (13) 

where we have have isolated the part that does not depend on the acceleration in the 
quantity 

J = 2N\\ (n • X) - hT 1 (X 2 ) (n • X") . (14) 

To make contact with the analysis in [19J for axisymmetric configurations, 
recall that an axisymmetric configuration can be described by cylindrical coordinates 
{R(t) , Z (t) , <f)} . The curve C is parametrized by the angle 0; the tangent vector is 
X' = <9X/<90, and the one- dimensional metric is h = R 2 . The velocity X is along a 



meridian, with N\\ = 0, and N = y R 2 + Z 2 . Moreover, in this case, n ■ X' = 0, n ■ X" = 
—N~ 1 RZ, and n ■ X = N X I 2 {ZR — RZ); it follows that the quantity J takes the simple 
form J = R~ l NZ, so that the mean curvature is K = R' 1 N-^ 2 [R(RZ - ZR) + N 2 Z}. 

We now cast the energy F[X] given by (J2J as an action decribing the dynamics in 
time t of a curve parametrized by u, X(i, u): 



F[X] = J dtL[X,X,X] , (15) 
where the Lagrangian functional is given by 



L[X, X, X] = j> du N \fh JLj (_ n -X+ jf + ^(-n-X + J 



a - -n-X 

3 



(16) 
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Note the dependence of the Lagrangian on time derivatives of the field variables X: it 
is quadratic in X, the acceleration of the curve, and it is not possible to eliminate this 
dependence using integration by parts within the action. It depends on the velocity X 
through its appearence in N, J and n. In all but the last term, the position X appears 
only through its derivatives with respect to the parameter u along the curve, X' and 
X"; note that there is an implicit dependence on X of N and n. 

The Hamiltonian formulation of the dynamical system defined by the action (JT5j) 
is somewhat involved. As an aid to establishing our bearings, it will be useful to first 
understand the Hamiltonian formulation of a simpler system, without higher derivatives 
to contend with. 

3. Soap bubble 

In this section we will provide a Hamiltonian formulation of the equilibrium of a soap 
bubble, in which surface tension is pitted against the internal pressure. Here the energy 
is simply proportional to the area of the bubble, subject to a constraint fixing the 
enclosed volume. The energy involves no higher than the first derivatives of X with 
respect to the parameter t. This system illustrates, at the Hamiltonian level, the 
consequences of the presence of a local symmetry - reparametrization invariance - 
namely the presence of constraints on the phase space variables. The reader should be 
warned, however, that the canonical structure is very different from the one presented 
below for a lipid vesicle. 
The energy is 

F[X] = a A - PV , (17) 

where a is the surface tension, and P the pressure excess inside the bubble. The 
equilibrium condition obtained from the vanishing of the first variation of F\X] is given 
by the classical Laplace- Young equation, of second order in derivatives of the shape 
functions, 

aK-P = 0. (18) 

Equilibrium configurations are surfaces of constant mean curvature. 

We consider the energy ()17j) as an action, decribing the motion in the parameter 
t of a curve parametrized by u, X(i, u). Decomposing the energy with respect to the 
coordinates {t, u}, we have 

F[X] = fdt L[X,X] , (19) 

where the Lagrangian functional is 

L[X,X] = jduNVh ^r-~n-Xj . (20) 

We emphasize that this Lagrangian depends at most on the velocities via iV and n 
related respectively to X by (JJJ) and (JTTj). The factor of \fh appearing in (|2U)l ensures 
invariance of L with respect to reparametrizations of the curve C. 
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Our first step towards a Hamiltonian formulation of the problem is the derivation 
of the canonical momentum p conjugate to X. A direct approach would be to define 
the canonical momentum p via the functional derivative p = 5L/5X. However, with an 
eye towards the higher order generalization for fluid vesicles, it is preferable to extract 
the canonical momentum from the first variation of the action 1)19)1; we do this in some 
detail. Under an infinitesimal deformation of the embedding functions X — ► X + SX, 
the first variation of the action ([19)1 . can be written as 

5F[X] = J dAE-5X + ^«p-5X, (21) 

where E denotes the Euler-Lagrange derivative. This expression allows us to read off 
the canonical momentum p from the second term. To carry out the first variation, we 
use the well known variational expressions (see e.g. |22j) 

5 A = JdAe a - d a 5X , (22) 

SV = - fdA [n • 5X + (n ■ X)e a • d a 5X - (e a ■ X)n • d a 8X] . (23) 
3 J 

Neither expression depends on the parametrization. For the first variation of the energy 
()19|) it follows that, integrating by parts it can be cast in the form (j21j) . 

5F = JdA (-<xV a e a - Pn)-<5X+ J dAV „ j<xe a • 5X - ^ [(n • X)e a - (e a • X)n] ■ <5X j .(24) 

Now we use the Gauss- Weingarten equations for the surface £ 

V a e 6 = - K ab n , (25) 
V a n =K ab g bc e C} (26) 

and Stokes' theorem as applied to a surface vector field V a which states that, for a 
region of the surface enclosed by the curve C, 

J dA V a V a = I duVhl a V a , (27) 

where l a is the normal to C on the surface S (note that 1 = l a e a ). Therefore we obtain 

5F = JdA (aK - P) n • 5X + j> duVh jal • 5X - ^ [(n ■ X)l - (1 ■ X)n] ■ 5x| . (28) 

Comparison of (|21jl and ()28j) identifies the Euler-Lagrange derivative as E = En, where 
E is given by the l.h.s. of (|18j): the Euler-Lagrange derivative is normal to the surface. 
For the canonical momentum we read off 



p = vfijri-£[(n.X)l-(l.X)n]J 



= a Vhl-^XxX' . (29) 

The momentum transforms as a density under reparametrizations of the curve C, as it 
should. It is normal to the curve. In this sense, it differs in an essential way from X 
which is not necessarily normal. The contribution to p from area is proportional through 
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cr to the densitized unit vector 1. The addition to p due to the volume constraint is 
independent of X, which introduces a component normal to the surface. 

Note that neither N nor N» are canonical variables for the soap bubble: they cannot 
be expressed as functionals of X and p. They should not therefore appear within the 
canonical framework. For the same reason, the normal n (unlike the tangent vectors t 
and 1) is also not a canonical variable, and has no place within this framework in this 
model. For fluid membranes, on the other hand, N, N» as well as n will be identified 
as functions of the canonical variables. 

We have now identified the phase space variables appropriate for a Hamiltonian 
description of the equilibrium configurations of a soap bubble: the position of a given 
curve X and its conjugate momentum p, as given by ([29)1 . The canonical Hamiltonian 
is obtained via a Legendre transformation of the Lagrangian (|20jl with respect to X as 



#o[p,X] = f dup-X-L[X,X] . (30) 

It is immediate to see that it vanishes, 

F o [p,X] = 0. (31) 

Its vanishing is a direct consequence of the presence of reparametrization invariance. At 
first this might appear to indicate some inconsistency. However, a local symmetry such 
as reparametrization invariance also implies constraints. These constraints will generate 
the evolution. Let us examine these issues one at a time. Reparametrization invariance 
renders it impossible to invert for the velocity in terms of the momentum. To see this 
explicitly one considers the Hessian 

5 2 L aVh 

HiA = : :— = HiUn . (32) 

Hij is degenerate: its determinant vanishes. This is a feature of any Hamiltonian system 
with local symmetries. It means that at any value of the parameter t the phase space 
variables are not all independent; there are constraints. We use the definition of the 
momentum (|29|) to identify the constraints as functions of the canonical variables: p is 
orthogonal to C, and the function of the canonical variables (p + P X x ~X' /3)/aVh is a 
unit vector. We thus find that, at each value of t, the phase space variables must satisfy 
the primary constraints: 

C, = j p - -XxX'l -a 2 h = (33) 




C 2 = p • X = . (34) 



The only freedom p possesses independent of X is a direction. 

The constraints possess a geometrical role in phase space as generators of 
infinitesimal reparametrizations. In particular, whereas C 2 generates reparametrizations 
of the curve C itself, C\ generates motions off the curve C onto the surface. This can 
be shown by considering the Poisson bracket of the constraints with the phase space 
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variables, where, for two arbitrary functions of the phase space variables / and g, the 
Poisson bracket is 

~5f Sg 



{f,9} = f 



du 



5X 5 P if ~ 9) 



(35) 




It should be mentioned that analogous constraints with P = occur in the Hamiltonian 
formulation of a relativistic bosonic string (see e.g. [T5j). 

The Hamiltonian that generates the motion of the curve C is proportional to the 
constraints themselves, 

H = I du (Aid + A 2 C 2 ) , (36) 

where Ai = \\{u 1 1) and A 2 = A 2 (u,t) are two arbitrary Lagrange multiplier functions 
that enforce the constraints. The time derivative of a phase space function / is given 
by its Poisson bracket with this Hamiltonian, 

/ = {/,#}• (37) 

In particular, it is important to note that the time derivative of the constraints (J3H|) 
and does not generate any new constraints. 

Hamilton's equations are given by (see the Appendices for useful formulas) 

X = ^ = 2Ai (p + ^X x X'^ + A 2 X' , (38) 

*--i-f^(- + ^*)-f[^ 

- (2a 2 \ 1 X' + \ 2 p)'. (39) 

These two couple first order equations reproduce the equilibrium condition E = En = 0, 
where E is given by the l.h.s. of the Laplace- Young equation (|18|). that follows from 
extremization of the energy. The first equation identifies both the momentum p and 
the Lagrange multipliers Ai,A 2 . Using the constraints, one finds that the Lagrange 
multipliers are proportional to the lapse and shift functions appearing in (|Bj): 

A ' = i^r X ' = N *- (40) 

These expressions for the Lagrange multipliers indicate that, given the curve C, the 
surface we generate depends only on the normal 1: the specific linear combination of 1 
and t entering the velocity X is gauge - it is not determined by the Hamilton equations. 
For the soap bubble, reparametrization invariance manifests itself in a freedom within 
the velocity X (for a fluid membrane this will become a freedom within the acceleration 
X). The definition (|29jl of the momentum is reproduced when we substitute PU|) for the 
Lagrange multipliers back into ([38)1 . Using this information in the second of Hamilton's 
equations (jSHJ), one finds the following relationship: 

-^£n = p+^|, (41) 
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where E is given by the left hand side of (|18|). To check this, one can use the geometrical 
expressions given in Appendix A. 

Hamilton's equations (jSHJ), ((HSj) , together with the constraints (jSSI), provide 
an alternative formulation of the single Euler-Lagrange equilibrium equation ()18|) . Once 
we have confirmed that this set of equations is equivalent to (fTH|) . we can forget about 
the latter. 

The basic strategy to construct an equilibrium configuration is to choose a closed 
curve in space defined by its position X(w, 0). We choose a momentum p(w, 0) consistent 
with the constraints (JHHj) and (|34"j). This corresponds to a direction. Now we let X, p 
evolve according to the Hamilton's equations (fHKj) and ffTTTTf) for any Ai and A2 of our 
choosing. An equilibrium configuration X(u, t) will be generated. 

In practice, when we discretize these equations, we will find that the canonical 
variables will stray off the constrained surface in phase space; the formalism itself 
suggests a strategy to nudge them back. One simply introduces explicit damping terms 
proportional to the numerical value of the constraints on the r.h.s. of Hamilton's 
equations. It is also possible to tune Ai and A 2 so as pace the evolution most 
appropriately. 

Unless we have been extremely lucky, we will find that the initial data we have 
prescribed generates a surface that does not close smoothly. This is not surprising: 
after all, the unique non-self-intersecting geometry is a sphere of radius i?o = 2<r/P. 
One needs to tune the initial conditions to produce this sphere. 



4. Phase space for the fluid membrane 

We turn now to the Hamiltonian formulation of the evolution of a curve described by 
the action ()15j) . We first identify the phase space variables. Because the Lagrangian 
(fTH|) depends on the acceleration X, in addition to the position of the curve X and its 
conjugate momentum p, the phase space will need to be enlarged to include the velocity 
X and its conjugate momentum P. (A pedagogical introduction to the Hamiltonian 
formulation of higher order systems is given in the Appendix of [19J.) The momenta are 
defined by the functional derivatives 

-£-*(£)■ 

Using the expression (|13|) for the mean curvature K, one immediately obtains for the 
momentum P conjugate to X, 

p = -^(nK + P)n. (44) 

P is normal to the surface, and transforms as a density under a reparametrization of the 
curve. The part due to bending is proportional to the mean curvature. We anticipate 
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that this expression will be the content of one of Hamilton's equations for this system 
(see (J?nj) below). Note that P depends on X, X, and X. 

The functional derivatives involved in the determination of the momentum p are 
somewhat less simple. A more convenient approach, generalizing that used in the 
previous section for the soap bubble, see (|2*Tjl . focuses on the first variation of the 
energy. For a system that depends on the acceleration, (}2~T]) generalizes to 

5F[X\ = J dA E • 5X + j> du (p • 5X + P • 5X) , (45) 

where E is the Euler-Lagrange derivative. Therefore, just as we did for a soap bubble, 
it is possible to read off the momenta from this expression. To carry out the variation, 
we use variational results described e.g. in (221, such as (1221). (1251). together with 



SK 



n • V 2 5X + 2K ab e a ■ d b 5X . 



(46) 



The Euler-Lagrange derivative is normal, E = En, where E is given by the l.h.s. of the 
shape equation (J3J). Let us focus our attention on the second term of (|4*HJl . Integrating 
by parts and using Stokes' theorem (|27jh one finds that 

{V a K)n + {-K 2 g ab - KK ab " 



jdu (p • 5X + P • SX 



= I duVhl a { 
+ (3{Kg ab - K ab 



e 6 



ie 6 



<re 



[(n-X)e a - (e a -X)n]} ■ 5X 



- j duVhl a (KK + p)n ■ V a 5X . 



(47) 



The terms in parenthesis spanning the first two lines contribute to the momentum p. 
The third line need some massaging to bring it into the form (j45j) : we use the definition 
of the velocity ©, to decompose l a V a SX = N~ 1 (SX - N\\dXf). The third line of (gTj) 
thus contributes to both P and p: 

- j du Vh(nKn + (3n)l a V a 5X = <fduP-(5X + N\\ (5X)') 

= ^ du [P ■ 5X - (N {l Py ■ 5X] . (48) 

The definition (|44|) of the momentum P is recovered, and we identify the momentum p 
conjugate to X as 

p = Vhl a 

+ aVh\ - ^rX x X' + (iViiP)' . (49) 



K 



-K 2 g° 
2 y 



KK ab ) e b + (V a K)n 



+ p(Kg ab - K ab )e h 



ab \ 



The momentum p also tranform as a density under reparametrization of the curve. Note 
that it depends on X, X, X and X with three dots. For k = (3 = the momentum for a 
soap bubble (J2~9"j) is recovered. We also anticipate that one of Hamilton's equations will 
reproduce this expression (see (j8()jl below). Geometrically, the momentum p given by 
(J49|) is very simple; p is also related in a very direct way to the conserved stress for fluid 
membranes derived by two of us in ^3]- To illustrate the relationship, it is convenient 
to set P to vanish. The conserved stress for the energy (0) is 



r 



KK° 



-K 2 g ab 
2 y 



e b - (V a K)n 



+ (3 (K ab - Kg ab ) e b 



ere" 



(50) 
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Its divergence gives the shape equation as a conservation law: 

V a P = £n = 0, (51) 

where E is given by the l.h.s. of (@J). We exploit the completeness of the surface tangent 
vectors g ab = l a l b + t a t b , where l a = 1 • e a , and t a — t ■ e a to decompose the stress as 

P = f l a + g t a . (52) 

The canonical momentum p is given by (minus) the projection of the stress along the 
normal to the curve l a , up to a total divergence: 

p = -y/hf + (N\\P)'. (53) 

Note that integrated over the closed curve: 

du p = - j> duVhl a f a . (54) 

The integrated momentum is the total force acting on the region enclosed by the curve; 
if there are no external forces acting and P = 0, this integral will vanish. 

One of Hamilton's equations will be an equation for p (see (181)) below). Modulo 
the remaining canonical equations, this equation should reproduce the conservation law. 
Using the decomposition of the stress tensor (|52j) in (J51|) . we find 

NVhVJ a = d t (Vhf) + (JVg - N\\Vhf)' . (55) 

where we also make use of the identities JHSl), ([100|l given in Appendix A. If we now use 
([53]). we find that the (densitized) shape equation can be cast as a continuity equation, 

y/gEn = yfgV a f a = -p + m' = , (56) 

with a current m given by 

m = «9 t (iV||P) +iVg + iV||p-iV||(iV||P) / . (57) 

If now we restore a non vanishing pressure P, the continuity equation ()56|) acquires a 
source term, see (jSIjl below. 

5. Hamiltonian and constraints 

The canonical Hamiltonian is given by a Legendre transformation of the Lagrangian 
L[X, X, X], as given by IjlfiJ). with respect to both the velocity X and the acceleration 
X, 

# [X,p;X,P] = jdu (P-X + p-X) -L[X,X,X], (58) 

where the definition of the momentum P is used to express the acceleration X in terms 
of the phase space variables. It should be emphasized that the term p ■ X is left alone; 
the only place p appears in Ho is in the first term. The only tricky bit is the elimination 
of X. To do this, first note that from the definition of P (I44j) . and the decomposition 
of the mean curvature ([130. we have 

P ■ X = ^(n ■ X) 2 - ^(U + /W 2 )(n ■ X) ; (59) 
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the quantity J, defined in (|14|). does not depend on the acceleration. Setting 
L[X, X, X] = I du £(X, X, X) , 

we find that 



P ■ X - £(X,X,X) 



VhN 



K 



fn-X) 



-J 







J 



a + -n.X 



(60) 



(61) 



2N Ay ' 2iV 4 " N 2 
the acceleration has been isolated in a single term. We can express this term with 
respect to the phase space variables by inverting the definition for (n • X) 2 : 



N 3 



«(n • X) - (kJ + I3N 2 



n 



(62) 



implies 



fn-X) 



N 6 
K 2 h 



P + 



~W Viv 2 



-, 2 



J + /3 



ii 



(63) 



Substituting this this expression into 1)61)1 . we identify the canonical Hamiltonian as 



# [X,p;X,P] = fdu [p-X + 



N 3 
2nVh 



P 2 +[J + N' 



P n + 



2k 



aNVh + -NVh(n ■ X)] . 
3 



(64) 



Unlike its soap bubble counterpart, this Hamiltonian does not vanish at the level of its 
definition, although it turns out that it does vanish anyway, see below. It is quadratic in 
the higher momentum P; it is linear in p — indeed, p only appears in the first term. The 
dependence on X is rather more complicated: linearly in the first term through p-X, and 
non-linearly everywhere else through its appearance in N, n and J. The dependence on 
X is similarly complicated; it appears always through its spatial derivatives with respect 
to the parameter u along the curve C with the exception of the last term proportional 
to the pressure P in which X appears directly. While this dependence does not break 
translation invariance, it signals that p will not be a divergence. 

The reparametrization invariance of the energy manifests itself, at a Hamiltonian 
level, in the presence of constraints: the phase space variables are not all independent at 
a fixed value of the parameter t. Primary constraints appear because it is impossible to 
invert for the higher momentum P in terms of the acceleration X, just as, in Sect. 4, it 
was impossible to invert for the momentum p in terms of the velocity X. Note, however, 
that this did not hinder the completion of the Legendre transformation. The primary 
constraints are identified by the null eigenvectors of the Hessian of the Lagrangian with 
respect to the higher derivatives 

£ 2 L[X,X,X] Vh 

1 = ~ K Jp n i n 3 ■ I 65 ) 

X and X'; we therefore identify the two primary 



n 



6X*6Xi 

There are two null eigenvectors 
constraints 

d = P • X = , 
c 2 = P ■ X' = . 



(66) 
(67) 
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Their geometrical content is simple: the momentum P is normal to the surface. 

The Hamiltonian function that generates the motion is given by adding the primary 
constraints to the canonical Hamiltonian 1)64)1 . giving the total Hamiltonian 

H[X, p; X, P] = H [X, p; X, P] + j> du {\ X C X + A 2 C 2 ) , (68) 

where Ai = Ai(w,t),A2 = \2{u,t) are arbitrary Lagrange multipliers that enforce the 
constraints. 

The Poisson bracket appropriate for this higher derivative system is, for any two 
functions on the phase space f,g, 



{f,9} = f 



du 



8f Sg 6f 6g ( 

5X ' 6P 8X ' 5p U ~ 9) 



(69) 



The time derivative of any function on phase space / is then given by its Poisson 
bracket with the total Hamiltonian 



f={f,H}. (70) 

We have identified the primary constraints on the phase space variables. For 
consistency, these constraints must be conserved in time. This requires us to introduce 
new secondary constraints. First, note that the primary constraints are in involution 
under the Poisson bracket, i.e. {C^C^ = Ci- The conservation in time of the 
primary constraints (Jfi7)j) and (|H7)l implies the existence of the secondary constraints, 
with H = § du Ho, 

5 1 = H = 0, (71) 

5 2 = p • X' + P • X' = . (72) 

A detailed proof that these constraints follow from the conservation in time of the 
primary constraints is provided in Appendix B. As we will show below, a far simpler 
derivation follows directly from reparametrization invariance. 

Among the secondary constriants is the vanishing of the canonical Hamiltonian, 
the hallmark of reparametrization invariance. Together, the two secondary constraints 
completely fix the two tangential components of the momentum p, p • X and p ■ X'. 
There are no other, tertiary constraints - the four constraints we have identified are in 
involution under the Poisson brackets. 

The secondary constraints are the generators of reparameterizations. In particular, 
5*2 generates reparametrizations tangential to the curve C, as follows from the fact that 
if we denote by Z any of the phase space variables, 

5Z = {Z, S 2 } = 71 . (73) 

On the other hand, S\ = —H,q generates reparametrizations off the curve C, as it will 
become clear in the next section, where we consider the evolution it determines. The 
details of the model will enter through Ho. 

Conversely, one can show that reparametrization invariance implies all four 
counstraints: the primary constraints ()66)) and (jBTj) as well as the secondary constraints 
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(J7TJ), ()72j). To see this, let us re-examine the expression for the first variation 
of the energy. Consider first a temporal reparametrization 5X = lu(u, t)X, where the 
parameter u is some function of the parameters w and t. In equilibrium, we have for 
the right hand side of 

y dw (p • SX + P • SX) = ^dwcj(p-X + P-X) + jduCoV -X, (74) 
On the other hand, for the left hand side of (|45p. we have 

£F[X] = jduuoL. (75) 
We conclude that 

y cj^o + uC x = . (76) 

Because u and u; are completely arbitrary, we conclude that the invariance of F under 
a redefinition of time implies the constraints ()66|) and (|71|1. In exactly the same way, 
reparametrization invariance tangential to the curve implies the constraints (jBTjl and 
(J72J1 . The stucture of the constraints depends only on the number of derivatives; it is 
independent of the details of the model. 

This analysis also clarifies the role of the primary constraints. In the transformation 
for the velocity induced by SX = o>X, SX. = uX + uX, the second term is a scaling 
of X. The constraint C\ is the generator of this scale transformation. An equal and 
opposite scaling of P is required to preserve its Poisson bracket with X. 

A result of reparametrization invariance is that the phase space variables provide 
a redundant description of our system. The four constraints do, however, permit us 
to identify the independent degrees of freedom: these are position of the curve X and 
its velocity X as well as the normal components of P and p (or, equivalently, X and 
X® (with three dots)). Whereas both X and X can be freely specified, the primary 
constraints kill the tangential components of P and the secondary constraints fix the 
tangential components of p in terms of the remaining canonical variables. 



6. Hamilton's equations 

In this section, we present the complete set of Hamilton's equations that follow from 
the Hamiltonian (|68|). We find it convenient to consider the equations in this order: 



5H ■ 5H 5H 5H 

9,X=-; d t X=-; ftP = -^; *P = - m - (77) 

There is one equation for each phase space variable; the equations are first order in time. 

(1) The first equation is 

d t X = 5H/5p = X . (78) 

Note that the only dependence on p in the Hamiltonian is through the linear term 
p • X appearing in the canonical Hamiltonian H (J64)) . Although it may appear to be a 
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mere identity, its content is the identification of the canonical variable X with the time 
derivative of X. 



(79) 



(2) The second of Hamilton's equation <9 t X = 5H/5P is (see Appendix B): 

d t X = X = ^=P + I J + N 2 - ) n + AiX + A 2 X' . 
KVh V K J 

Note that the r.h.s. does not depend on the momentum p conjugate to X. The content 
of this equation is to identify the momentum P in terms of X, X and X, reproducing 
the expression (J44j) . This equation also identifies the form of the Lagrange multipliers 
Ai and A 2 . As we show in the next section, they are proportional to the tangential 
components of the acceleration X. 

(3) For the third of Hamilton's equation d t P = —SH/S'X., we obtain (see Appendix C): 



t P= - p + 2/T 1 (p ■ X')X' + 2(JV||P)' + {h^N- 1 
+ 2/T 1 (P • X") X 



2NMI ■ X') - X 2 (1- X") 



2kVK 



K 



2k 



1 - -X x X' . 

3 



Ai} P 

0) 



This equation identifies the momentum p in terms of X, X, P and P. It is independent 
of A 2 . This presentation of p is very different from the form written down in (J4"9~|) where 
it is cast as a functional of X, X, X and X with three dots. A little work is required to 
establish that indeed the two are the same p. We will return to this issue in the next 
section. 

At this point, it is worth emphasizing that the secondary constraints identified in 
the previous section, (|7Tj) and (|72^) . specify the tangential part of p in terms of the 
remaining phase space variables, and initial data will need to be chosen accordingly. 



rri 



(4) Finally, the fourth Hamilton equation dtp = —5H/ 5X is 
d t p = p 

where the current m is given by 

A^iiP 2 + — A^V||(P • n) - 2Nh- l (P ■ X') + VhN» 



-NVhn 

3 



(81) 



m 



3A^ 
2nVh 



K 



2k 



a 



+ 



+ 



h-" 2 — P 2 
2k 



h~ 2 ± 2 (P ■ X") + 2/r 1 AT||(P • X') - h- 1/2 N (— 



a 



X' 



Ny/h 



h-^in x X ■ X") - 2JVnfn xtx'l 



- A 2 P + -Xx X 

3 



12) 



Modulo the other Hamilton's equations, which establish the relationships between the 
dynamical variables appearing on the r.h.s., this is the shape equation (J3J) cast in 
canonical form. We will show this explicitly in the next section. Note that (|81|) 
assumes the form of a continuity equation. The origin of the source term is the explicit 



Hamilton's equations for a fluid membrane 



18 



dependence of the Hamiltonian if on X appearing in the volume when the pressure P 
is non-vanishing. 

To summarize, Hamiltons equations (fTHJ) . (fTH|) . (|KU|) and (JEIJ) together with the 
primary constraints (Jo^j). (foT)). and the secondary constraints (fTTj) . (|72jl. provide an 
alternative formulation of the shape equation Pj). Just as it does for the Einstein 
equations in General Relativity, this formulation calls for a change of perspective. We 
started with with a surface satisfying the shape equation. There is now no surface, only 
a curve. 

The recipe for reconstructing a surface is simple: (i) We choose any closed curve in 
space (X(0,w)). (ii) We are free to specify its velocity X(0, u). At this point we have 
defined a basis in space adapted to the curve given by the tangent to the curve, X'(0, u), 
X(0, u) and their cross product, proportional to the unit normal n(0, u). The remaining 
initial data on this curve must be consistent with the constraints. (Hi) To satisfy the 
primary constraints ()66|) and (jpTf) . the momentum P is chosen to be orthogonal both 
to X and to the tangent to the curve X'. (iv) The momentum p is constrained by the 
secondary constraints (fTTj) and (|T2j) which determine its two tangential components in 
terms of the initial data we have already specified. (Note that p must also be consistent 
with the integrability condition associated with a closed curve given by (jS3)). (v) We 
choose appropriate values for the two functions Ai and A 2 — this corresponds to a choice 
of coordinates. As we show in the next section, they are proportional to the tangential 
components of X. (v i) We now let the system evolve according to Hamilton's equations. 
An equilibrium surface X(w, t) satisfying the shape equation will be generated. 



7. Shape equation from Hamilton's equations with constraints 

In this section, we show how the Hamilton equations, together with the constraints, 
combine to reproduce the shape equation (jlj). 

The second Hamilton equation (fTHj) identifies P and the form of the Lagrange 
multipliers. We dot (|T§j) with the normal n, and we use the primary constraints 
(jnij), to obtain 

N 3 P ■ n = K y/h(X- n- J) - N 2 Vh(3 . (83) 

As the primary constraints tell us that P is purely normal to the surface, using the 
expression (|T3*j) for the mean curvature K, we recover the expression (jUj) for P. 

To obtain the form of the Lagrange multipliers we dot (|79|) with X and X', and use 
again the primary constraints (jMj) . (joTf) . so that we find 

X-X , . 

^ = ~%2~ ' ( 84 ) 

AT 

A 2 = fe~ 1 X / -X-JVnAi= nxX-X. (85) 

11 v^X 2 1 ; 

The first tells us that Ai has the geometrical meaning of the affine connection along the 

integral curve of X; the second that A 2 is proportional to the tangential component of 
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the acceleration X orthogonal to the velocity X. Therefore, for this higher derivative 
system, the components of the acceleration tangential to the surface are arbitrary; they 
are the gauge part of the evolution, whereas the normal component is proportional to the 
momentum P. In the case of a soap bubble, the Lagrange multipliers were proportional 
to the components of the velocity - this is its natural higher order generalization. 

The third Hamilton equation (}80|) identifies the form of the momentum p, as given 
by (|49j). To see this is not immediate. To facilitate the comparison, we consider first 
(JH3). We use the Gauss- Weingarten equation (J2SJ), and the expression (jUj) for P, to 
obtain for p 



P 



\[hl a 



-K 2 g ab 
2 y 



2KK ab ) + (3{Kg ab - 2K 



ab 



e 6 - P + 2(iV||P)' 



P + aVh 1 - -X x X' . 

3 



Now using the completeness relation on the surface g ab = l a l b + t a t b , we have that 

1 



K ab l a e b = Kl 



hN 



(n ■ X")X - (n • X')X' 



(86) 



(87) 



and inserting this expression in (|8fi)l we obtain the alternative expression for p 



p = sfh K 1--K 2 1 - 2h- 1 N- 1 K 



n-X")X-h- 1 N- 1 (n-X.')X' 



- (3 {VhK 1 + 2h~ 1/2 N~ 1 [(n • X")X - (n • X')X'] } - P + 2(JVj|P)' 



" Vh a \N / 



P + aVh 1 - -X x X' . 

3 



Using the information gathered from the second Hamilton equation (|79)l. i.e. (144)1 for 
the momentum P, and (|54^1 for the Lagrange multiplier Ai, it is a simple matter of 
algebra to see that the third Hamilton equation (jHOj) coincides with the expression (JHHJ) 
for p. 

Finally, the fourth Hamilton equation (|5T|l is (|3B^I. that we have derived in Sect. 4 
from the conservation of the stress tensor at equilibrium. To see this, requires using the 
specific form P,Ai,A2 in the expression ()82)1 for the current m, and comparison with 
the right hand side of the alternative expression (|57)). 



8. Concluding remarks 

In this paper, we have presented a Hamiltonian formulation of the shape equation 
for arbitrary, not necessarily axisymmetric, configurations of lipid vesicles. In this 
description the surface geometry is reconstructed from a closed curve. This is possible 
because the conserved internal stresses which underpin the equilibrium geometry are 
transmitted across closed curves; the shape equation is a conservation law for these 
stresses. 

This formulation of the shape equation offers a new approach for constructing 
equilibrium configurations based on the numerical solution of Hamilton's equations. 
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Whether its implementation turns out to be tractable, in practice, remains to be seen. 
There is no problem generating equilibrium surfaces: any initial data consistent with the 
constraints will do this. The subtlety is identifying initial data that are consistent with 
a surface which closes smoothly. In general, the curve will self-intersect and singularities 
will arise. The initial data will thus need to be tuned accordingly. This problem is a 
global one. Progress is likely to be limited until a second problem is tackled: what is 
the most appropriate choice of Lagrange multipliers for any given initial data? This 
will fix the parametrization of the surface. There may not be a single choice which is 
appropriate everywhere. 
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APPENDIX A 

In this Appendix we collect various geometrical formulae that are useful at 
intermediate stages of the calculations that we have performed. The Frenet-Serret 
equations for the curve C, seen as living on S are, with V s = t a V a : 

V s t a = kl a , (89) 
V s Z a = - k t a , (90) 

where k is the geodesic curvature for the curve. The Frenet-Serret equations for the 
curve C seen as living in space are 

V s t = fcil+ k 2 n (91) 
VJ = -ht + hn (92) 
V s n = -k 2 t-k 3 l (93) 

where ki — k is the geodesic curvature, and k 2 , k 3 are given by the projections of the 
surface extrinsic curvature by k 2 = —K a bt a t b , k 3 = —K a bt a l b . 

The evolution in the parameter t of the basis {X', 1, n} can be derived from the fact 
that partial derivatives commute, so that X' = (X)'. It follows that 

X' = {d u N + N\\k) 1 + (d u N\\ - h^Nk) X' + (n • X') n , (94) 

1 = — h" 1 (N' + JV|| fc) X' + iV -1 [(n • X) — iV|| (n • X')] n , (95) 

ri = — /i -1 (n • X') X' — A _1 [(n • X) — iV||(n • X')] 1 . (96) 

Note that 

d t Vh = -VhNk + Vhd u N l{ . (97) 
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The acceleration is 

X= (N ll -h' 1 Nd u N + N ll d u N ll -2h- 1 NN ll k) X' + (N+N {l d u N+N^k) l-(n-X) n.(98) 
We also note the identities 

VX = - k (99) 
y a t a = N~ 1 V S N. (100) 
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In this Appendix, we show how the secondary constraints (fTTj) . ()72|) follow from 
the conservation in time of the primary constraints (jHEJ), (|H7j) . The first thing to note 
is that the primary constraints are in involution under the Poisson bracket (|69jh 

{C 1 ,C 2 } = C 2 nO, (101) 

where the symbol ~ denotes equality modulo the constraints themselves. 
The time derivative of the first primary constraint is 



C 1 = {C 1 ,H}^{C 1 ,H } = P- 



X 



SHn 



(102) 



5P 5X 

We need to evaluate the functional derivatives of the canonical Hamiltonian functional 
(164)1 . For the functional derivative of the canonical Hamiltonian with respect to P, we 
have immediately that 

K 



11 



(103) 



In order to evaluate the functional derivative of H with respect to the velocity X, we 
hold X fixed, and consider the variation with respect to X of the geometric quantities 
that appear in the canonical Hamiltonian: 

5N =1-SX, 

5N\\ = h- 1 X'-SX, 

5n = -N- 1 (n-5X) 1. 

For the part of the mean curvature that does not depend on the acceleration X, defined 
as J in ([14)1 . it follows that 



N~ 



5 J = h- 1 {2(n ■ X')X' 

+ 2A|| n ■ (6X)' . 
Using these expressions we have 

2N(3. ri 



X 2 (l • X") - 2N\ { h(\ ■ X')] n - 2(n • X")x} • SX 



5X 



P 



3N 2 
2kVK 



p2 + 



K 



n) + - - a \fh 

2k 



l + 2h _1 (n-X')(P-n) X 



2[iV||(P • n)n]' + h^N-^P ■ n) [x 2 (1 • X" - 2N\\h(l ■ X)) 



n 



- 2h~ 1 (n ■ X")(P ■ n) X + -X x X 



(104) 
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We can simplify this expression using the completeness relation 5 U 
hr x X l X\ so that, modulo the constraints, (P ■ n)n w P. Hence 

SH 



22 



5X 



P 



3JV V + !^ ( P. n ) + <^_^ 



2k 



1 



-l j\r-l 



X 2 (l-X"-2A,/i(l-X)) 



P - 2/rVP ■ X") X 



2/T 1 (P-X')X' 
P 



2(iV||P)' 
XxX'. (105) 



and one finds that the conservation in time of C\ gives the vanishing of the canonical 
Hamiltonian density TCo, defined by H = § du TLq\ 

C 1 « -H = . (106) 

To see that this expression holds it is essential to recognize that 

J(P • n) « iV||(P • X') - (N\{P)' ■ X - h- l X_ 2 (P ■ X") , (107) 

where again one uses the fact that (P • n)n « P, and (AiP)' • X' w — Ah(P • X') . 

Similarly, the conservation in time of the second primary constraint (JBTj) gives the 
secondary constraint ()72j) : 



C 2 = {C 2 ,iJ}^{C 2 ,^o} = -X'- 



<5X 



X P' 



-s 2 = o 



(108) 
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In this Appendix we fill in the details of the derivations of the third and fourth 
Hamilton equations (jSUj) and ([81)1 . Consider first the third Hamilton equation, P = 
—5H/5X. We use ()105j) from Appendix B, and we obtain immediately 

SN 2 „, 



P 



,2/tv^ K 
+ /^iV" 1 [2iV||/i(l ■ X') - X 2 (l ■ X") 
P 

XxX'-AiP. 



2n 



l + 2(iViiP)'-2/r 1 (P-X)X 



P + 2h~ 1 (P -X")X 



(109) 



Using the secondary constraint ()72)1 in the last term of the first line gives the third 
Hamilton equation in the form given in the text, (|80|) . 

Let us turn to the derivation of the fourth Hamilton equation p = —5H/5X. Now 
we hold X fixed, and consider a variation of of the geometric quantities that appear in 
the canonical Hamiltonian w.r.t. X: 

sx' = (sxy , 

51 = - A" 1 A,, [n • (5X)'} n - h~ x [1 ■ (5X)'} X , 

Sn = (n~ x N\{1 - h~ x X') [n ■ (SX)') = N- x hr 1/2 (n x X)n • ($X)' , 

SN = - An 1 • (5X)' , 



SNn = h~ 



Al ■ (5X)' - Ay X ■ (5X)' 
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For the quantity J this implies 

5 J = 2Nh-\n ■ X')l • (5X)' + 2 [/r 2 X 2 (n ■ X") - iVy/T^n ■ X')] X' ■ (<JX)' 

n ■ (5Xy 
(110) 

Using these expressions, and exploiting again the fact that (P • n)n ^ P, we find that 
the fourth Hamilton equation is 



N» 



2_i(l . X') - AA- 1 AT | |/ i - 1 X 2 (l ■ X") - 2iV||/r 1 (X / ■ X') + /T 2 X 2 (X' ■ X") 



+ 

- h^X 2 n • (<5X)" 



d u 



3N 2 



2(3 



N\\P 2 + -t-NN\\CP ■ n) - 2Nhr l (P ■ X') + VhNn — - 



K 



2k 



+ 



+ 



/T 3/2 ^P 2 _ /,-2x2( P . X ") + 2/>- 1 iV||(P ■ X) - h- l ' 2 N - a 



a 



u 



1 



iV 



/i _1 iV||X 2 (l • X") - 2A|(1 • X') + 2iViV||/r 1 (X' • X') - /r 2 iVX 2 (X' ■ X") 
- (x 2 /i- 1 P)'-A 2 P + ^Xxx|-^iVv / /In. 
This gives the fourth Hamilton equation in the form given in the text, ()81j) . 



P 

(111) 



Hamilton's equations for a fluid membrane 



24 



References 

[1] Lipowsky R, and Sackmann E 1995 (eds.) Handbook In Biological Physics vols 1,2 (Amsterdam: 
Elsevier) 

[2] Boal D 2002 Mechanics of the Cell (Cambridge: Cambridge U. Press) 
[3] Canham P 1970 J. Theor. Biol. 26 61 
[4] Hclfrich W 1973 Z. Naturforsch. C28 693 
[5] Evans E 1974 Biophys. J. 14 923 

[6] Willmore T J 1982 Total Curvature in Riemannian Geometry (Chichester: Ellis Horwood) 
[7] Svctina S and Zeks B 1989 Eur. Biophys. J. 17 101 
[8] Seifert U 1997 Adv. in Phys. 46 13 

[9] Bozic B, Svetina S, Zeks B and Waugh R 1992 Biophys. J. 61 963 
[10] Wiese W, Harbich W and Hclfrich W 1992 J. Phys.: Condensed Matter A 1647 
[11] Miao L, Seifert U, Wortis M and Dobcreincr H G 1994 Phys. Rev. E 43 5389 
[12] Ou-Yang Z C and Hclfrich W 1987 Phys. Rev. Lett. 59 2486 
[13] Ou-Yang Z C and Hclfrich W 1989 Phys. Rev. A 39 5280 
[14] Capovilla R and Guvcn J 2002 J. Phys. A: Math, and Gen. 35 6233 
[15] Deuling H J and Hclfrich W 1976 J. Physique (France) 37 1335 
[16] Gompper G and Kroll D M 1997 J. Phys: Condensed Matter 42 8795 
[17] Bowick M and Travcssct A 2001 Phys. Reports 344 255 
[18] Capovilla R, Guvcn J and Rojas E 2004 Class. Quant. Grav. 

[19] Capovilla R, Guven J and Rojas E 2005 Hamilton's equations for a fluid membrane: axial symmetry 
preprint 

[20] Miiller MM, Deserno M and Guvcn J 2005 Euro. Phys. Lett. 69 482 

[21] Guvcn J 2004 J. Phys. A: Math, and Gen. 37 L313 

[22] Capovilla R and Guvcn J 2004 J. Phys. A: Math, and Gen. 37 5983 



